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In the framework of quark-hadron duality, the cross section values of charmonium 
and doubly charmed baryon production in e + e~-annihilations have been estimated 
at interaction energy yfs = 10.6 GeV. It has been shown that the approach under 
discussion allows to resolve some contradictions between the experimental data and 
pQCD predictions. 



I. INTRODUCTION 

Until recent times, the calculation of J/ip production cross section was done within the 
so-called 5-approximation, which is based on the factorization theorem of QCD. In the 
approach under discussion, relative momenta of valence quarks of J/^-meson are not taken 
into account in the hard part of amplitude. One supposes that both valence quarks are on- 
shell and each of them carries a half of the total J/^-meson momentum. Nevertheless, the 
recent experimental results of the BELLE and BABAR Collaborations on J/^-production 
in e + e~ annihilations demand an essential revision of the calculation techniques based on 
QCD factorization theorem. Indeed, the cross section value of charmonium pair production 
in e + e~ annihilations estimated in the framework o 
the experimental data by an order of magnitude 
that such underestimation is caused by large fixed virtualites of the intermediate quark and 
gluon q 2 ~ ^- (Q 2 is a virtuality of the initial photon, see Fig. HJ, which occur within the 
^-approximation technique. 

Taking into account the relative motion of valence quarks allows to remove contradictions 



' the (^-approximation ll|| underestimates 
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% y|. The detailed analysis has shown 
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between theory and experiment, at least in principle. Calculations of the pair quarkonium 



4] 



are in a qual- 



production in e + e~ annihilations, obtained within a light cone formalism 
itative agreement with the experimental data of the BELLE and BABAR Collaborations. 
In these works, the 5-shaped wave function of quarkonium was replaced by a function f(x), 
which is "spread" on x, where a; is a momentum fraction of the quarkonium carried by a 
valence quark in the infinite momentum frame. Thus, the longitudinal "internal motion" of 
the valence quarks is taken into account 1 . This spreading leads to a decrease of effective 
virtuality of the intermediate quark and gluon, and, therefore, to an increase of the predicted 
cross section value. 

In spite of certain success in the description of charmonium pair production, uncertainties 
of the concrete model realization remain. These are uncertainties of wave function value, 
final c-quark virtualities, and standart uncertainties of running constant value a s , c-quark 
mass value m c , etc. 

The cross section of inclusive charmonium production in e + e~ annihilations (e + e~ — > 

nn 

charmonium + c + c) depends on model uncertainties to a lesser extend |a, l|, because the 
largest cross section fraction for such process is determined by the smallest virtualities of the 



intermediate quark and gluon (Fig. EL Taking into account the relative motion of quarks 
changes the the result less than 10% g], and a general dependence on choice of m c and a s 
values remains. 

In this work we perform a new analysis of four charmed quark production in e + e~ anni- 
hilations within the quark-hadron duality. 

II. EXCLUSIVE AND INCLUSIVE CHARMONIUM PRODUCTION 

It is reasonable to describe the cross section of four c-quark production in the process 

e + e~ — > cccc, (1) 

in the framework of perturbative QCD. A value of m c = 1.25 GeV can be chosen for the 
c-quark mass. This value is the same as that adopted for QCD sum rules. The strong 



1 It is worth to note that the problem of internal motion of valence quarks inside J/ V'-meson was discussed 
a rather long time ago. See, for example, works 0, where it was suggested to use the structure function 
of J/ip in form F c = x 22 (l - x) 2A5 . 



3 



coupling constant at the scale of m c is put to a s = 0.24. It is worth to mention that the 
choice of relatively small value of c-quark mass makes a use of quark-hadron duality more 
convincing, because the masses of observed doubly heavy hadrons lie approximately in the 
center of duality interval at this m c value. 

The cross section value of process JTJ) can be considered as a sum of cross section values 
of all process with four c-quarks. Leading order calculations with the chosen values of a s 
and m c at y/s = 10.6 GeV give the following cross section value of the process 

a{e + e- cccc) = 372 fb. (2) 

The cross section of pure electro-dynamical production of four c-quarks is only 1 2% of 
the above value. 

There are uncertainties in prediction for the process (JTJ) due to the choice of m c and a s 
values. The m c and a s used here lead to the most optimistic prediction for the cross section 
value of four c-quark production. This value can be considered as an upper limit of four 
c-quark yield. The performed estimation ^ is in accordance with the cross section value at 
Z°-boson peak 

n 

As it was noted earlier |1CJ| , for the processes 

e + e" -> J/ip + c + c, (3) 

e + e _ — > T] c + c + c, (4) 

the total cross section values estimated within the (^-approximation [f3,0] agree with the cross 
section value of production of color singlet cc-pair calculated for a some duality interval. 

Using values of wave function at origin |_R S (0)| 2 = 0.67 GeV 3 and the charmonium mass 
value 3 GeV, the following results have been obtained for J/ip and r\ c production at y/s = 
10.6 GeV: 

a(e + e~ -> J/ij + c + c) = 104 fb, (5) 
a(e + e- ^ r/ c + c + c) =40 fb. (6) 

Therefore, 

a{J/i>)/a{ Vc ) = 2.6. (7) 

A contribution of higher S- and P-wave charmonium excitations into the inclusive pro- 
duction cross section is about 50% Q|, consequently the total cross section value can be 
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estimated 2 as 216 fb. This value can be compared to the cross section value obtained in 
the appropriate duality interval of mass spectrum of color singlet cc-pairs produced in the 
process JH): 

2m c < m c5 < IrriD + A, (8) 

where m c = 1.25 GeV, A = 0.5 GeV, and mo is the experimental D-meson mass. Thus, 
one can obtain: 



2m D +A ^ (e+e - -> (cc) sing i ct + c + c) 

am. 



= 280fb. (9) 

a=o.5 GeV 



/2m c rf^cc 

This value should be compared with 216 fb. Therefore, for the case of inclusive charmonium 
production (JHJ and (J3J) for some A a fair agreement between the results of ^-approximation 
and the quark-hadron duality hypothesis can be achieved. For the color singlet pairs with 
a definite spin value (S = or S = 1) one can obtain: 



2m D+ A da(e + e- -> (cc)| nglet + c + c 



2m c 



dm c5 

am, 



' cc 



5=1 



= 204fb, (10) 

A=o.5GeV 



2m D+ A da ( e + e - (cc)? ng lct + C + c) * " 

dm 



70 lb. (11) 

/2m c dm c5 A=o.sGeV 

If one supposes that the spin-triplet cc-pair (S = 1) most probably transforms into a J /if) 
meson and the spin singlet cc-pair (S = 0) most probably transforms into r] c , then 

a(J/if))/a( Vc ) = 2.7. (12) 

This ratio is in a good agreement with one calculated in the frame work of the 8- 
approximation. 

Thus for the case of inclusive charmonium production, the both approaches under dis- 
cussion agree. Contrary, for a case of exclusive charmonium pair production the cross 
section value calculated within the quark-hadron duality hypothesis is larger than the 5- 
approximation prediction by one order of magnitude. Such a difference occurs because the 
quark-hadron duality allows to take into account the internal motion of quarks inside char- 
monium. Indeed, a quark and an antiquark produced in duality interval (jHJ and fused into 
charmonium according to the quark-hadron duality hypothesis can have different velocities. 



Here and below the process of one-photon annihilation is discussed. Two-photon annihilation e + e 
— > J I if) + c + c contributes less than 10% into the total cross section value (ll] |. 
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This decreases the effective virtuality of the intermediate gluon and quark, and increases 
consequently the predicted cross section value. 

The exclusive charmonium pair production cross section can be estimated by picking 
out the both cc-pairs produced in the process (jlj) in appropriate duality intervals and in 
color-singlet state. In our calculations, we will also distinguish spin states of the pairs (Si 
and S2, correspondingly). The calculation results are as follows: 

r 2 mD +A d 2 a (e+e~ -> (cc) J nglct + (cc) s t nglct ) 



2m c dm^ dm C c 2 



-dm c5l dm c52 



2m D +A d 2 a (e+e -> (cc)f^ glet + (cc) s t nglct 

-Gtm c5l dm c52 

2m c dm r r-, dm. 



L CC\ " "CC2 



ff 2m D +A d 2 a (e+e -> (cc)J nglet + (cc) s t nglet 

; ; -am C c 1 dm C c 2 

am r r, dm, 



J '2m, c 



Si =0,S2=0 

= 1.2fb, (13) 

A=o.5GeV 

Si=0,5 2 =l 

= 19.4fb, (14) 

A=o.5GeV 

Si=l,S 2 =l 

= 19.8fb. (15) 

" CC ^"" CC2 A=o. 5 GeV 
One can suppose that if Si = and S 2 = 1 (or, S% = 1 and S 2 = 0), then J/ip and r\ c are 

produced in most cases; if S\ — 1 and 5*2 = 1, then and P-wave state of quarkonium 

are most probably produced. Then the cross section of process 

e + e~ — > J/ip + charmonium (16) 

at y/s = 10.6 GeV can be estimated to be about 40 fb. This value can be compared with 
the experimental data on the charmonium pair production: 

<r(e + e - — > J/ip + charmonium) = 55 ± 10 fb (BELLE Collaboration), (17) 

cr(e + e~ — > J/tp + charmonium) = 44.3 ± 9fb (BaBar Collaboration). (18) 

One can conclude that the quark-hadron duality prediction does not contradict the ex- 
perimental data within experimental errors and uncertainties in m c and a s values. There is 
certaintly no reason to talk about the discrepancy of order of magnitude, as it is a case for 
the predictions of the ^-approximation jlj]. It is worth to note, that the calculation results 
(fT4^l and (fl~5|) indicate a large yield of P-wave states, which is comparable with the yield 
of rj c and r( c . It is also in agreement with the experimental data of the BELLE and BaBar 
Collaborations. 

In the work jl^] the quark-hadron duality hypothesis was combined with the Colour 
Evaporation Model predincting a yield of charmonium state with the total angular momen- 
tum J is proportional to 2 J + 1. Obviously, this model cannot be applied to the process 
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e + e~ — > J/ip + X because the quantum numbers of J /if) and X channel are strongly corre- 
lated. 

The above estimations allow to predict the cross section of process e + e~ — > J/ip + DD. 
Subtracting the cross section value of exclusive pair production process from the cross section 
value of inclusive charmonium production process one can obtain: 

a(e + e- -»• J/ip + D + D)~ 240 fb. (19) 

That value can be compared with that obtained by the BELLE Collaboration (lsj : 

a(e + e- -> J/ tp + c + c) = 870 ±26 fb. (20) 

Thus, our predictions underestimates the experimental results by several times. Moreover, 
within the model under discussion the cross section value for the process e + e _ — > J/ip + DD 
cannot be larger than the cross section value of four charmed quark production, i.e. it can be 
not large then 330 fb, and the essential fraction of this value should include doubly charmed 
baryon contribution. It is discussed in the next chapter. 



III. DOUBLY CHARMED BARYON PRODUCTION 

Let us estimate how often a cc-pair is produced in the color antitriplet state with a 
small invariant mass in a system of four c-quarks. Such an object can be considered as 
a candidate to create a doubly charmed baryon (see also jjllj). In the duality interval 
2m c < m cc < 2m D + A for a s = 0.24, m c = 1.25 GeV, and A = 0.5 GeV, one obtains the 
following cross section value: 

a (e + e~ -> (cc)- 3 /(cc) 3 + X) = 170 fb. (21) 

For the production of both cc- and cc-pairs in the above duality interval, the cross section 
value is 

a (e + e- -> (cc)g + (cc) 3 ) = 32 fb. (22) 

We assume that in the kinematical region under discussion a doubly charmed baryon 
and an antibaryon are created with high probability. The rest of the cross section value 
cr(e + e~ — * (00)3/(00)3 + X) (about 140 fb) should be considered as an upper limit for the 
cross section of the processes like e + e~ — > E cc + A c + D. 
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Let us discuss here a relation between the inclusive charmonium production and the pro- 
duction of doubly charmed baryons. The c-quark and c-quark produced in the intermediate 
gluon splitting have rather often almost the same velocities and this provides the minimum 
virtuality of intermediate gluon (~ 4m;?). Let us suppose that the color singlet cc-pair is 
produced with a small invariant mass. Either c-quark or c-quark in this pair is produced by 
the splitting gluon. For the sake of simplicity, let us suppose that it is a c-quark. Then the 
remaining c-quark (which does belong to the pair under consideration) is also produced by 
the splitting gluon. It means that there is a large probability for this c-quark to have a ve- 
locity close to the velocity of the color singlet cc-pair. This is why the remaining c-quark can 
create with c-quark from cc-pair an antitriplet state with a small invariant mass. Therefore 
such a ccc-system can create a charmonium, as well as a doubly charmed baryon. Within 
the ^-approximation, the relative inclusive yield of charmonia and doubly charmed baryons 
are determined by the ratio of squared wave functions at origin for charmonium and doubly 
charmed diquark, correspondingly. Unfortunately, the quark-hadron duality technique does 
not allow to definitely separate these yields for the case of inclusive production. 

Contrary, the exclusive process of charmonuim pair production and the process of pair 
production of doubly charmed baryons are well separated kinematically, and the problem 
under discussion does not arise in this case. 

IV. PROCESS e+e" -> J/ijjgg 
When the charmonia production in e + e~-annihilation is discussed, the process 

e+g- ^ J/^gg (23) 

cannot be ignored (see Fig. [2j). The detailed study of the process is worthy of a separate 
article, because the situation here is quite similar to one with the exclusive pair production of 
charmonia (JT6l) . The treating of the process (|T6l) within the ^-approximation leads to large 
fixed virtualities of intermediate quarks, as in the case of process (jTfil) . This means, that 
the (^-approximation leads to the underestimated value of the cross section of the process 
(jH)) , as well as (|T6|) . The value of effective virtuality predicted within the quark-hadron 
duality is less then the value predicted within the ^-approximation, and it leads to the cross 
section value increase. In this work we do not intend to study the process e + e _ — > J/ipgg 
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in details and only give the cross section value predicted in the framework of hypothesis of 
the quark-hadron duality. For A = 0.5 GeV and for the invariant mass of two final gluons 
larger than two masses of 7r-meson (m gg > 2M^), one can obtain: 

2mD +A da(e + e- -> (cc)f ; lct + gg) 

dm 



2m c dm, 



cc 



2m D+ A da(e+e -> (cc)* nglet + gg) 



= 1.96pb, (24) 

A=0.5GeV,m oo >2M^ 



dm r 



ij>j 
s=o 



= 1.00pb. (25) 

A=0.5GeV,m 99 >2M» 



hm c dm C c 
If one suppose that the pair with a spin value of S = 1 most probably creates J/ip-state, 
then the results ([23$ should be considered as the cross section value estimation for the process 
e + e _ — > J/ipgg. Of cause, this estimation is quite rough. In addition, our calculations show 
that the prediction value depends quite strongly on A. For example, at A = 0.1 GeV the 
cross section value of the e + e~ — > J/ipgg is about 0.6 pb, which is in agreement with the 
experimental data. To achieve a better reliability of the predictions, one can introduce an 
effective mass of gluon into calculations, as it was done in the works |lJ, LL5J. Therefore, the 
quark-hadron duality hypothesis does not contradict the experimental data. Contrary, the 
predictions of the ^-approximation 13| 3 are too small to describe the data. 



V. CONCLUSIONS 

Until recently, it was believed that the experimental data on J/ip production in e + e~- 
annihilation cannot be described within the perturbative QCD. This conclusion was based on 
calculation which does not take into account the motion of the valence quark inside the J/i/j- 
meson. Taking into account the "internal motion" by means of light cone wave function, 
as it is done in works 0|, or by means of quark-hadron duality, as it is done in the present 
work, allows to resolve some contradictions between the pQCD and the experimental data. 
One can certainly claim that the experimental data on the charmonium pair production 
in e + e~-annihilation, as well as the experimental data on charmonium production in the 
process e + e~ — > J/xjjgg, do not contradict pQCD predictions. 

Nevertheless, the experimental value of the inclusive cross section of J/ip production in 
the process e + e~ — > J/ip + D + D is still underestimates by theoretical predictions. Even our 

3 In the work the factor of 1/2 has been missed in the cross section calculation for the process e + e _ — > 
J/tpgg. The revised prediction value is about 0.35 pb instead of 0.7 pb. 
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very optimistic estimation of the cross section value (330 fb at 10.6 GeV) is approximately 
three times less than the experimental result. Moreover, we think that the value of about 
200 fb seems to be the more reasonable estimation of the process cross section. 

In conclusion it is worth to mention one more time, that the calculations done in this 
work within the quark-hadron duality hypothesis suggest, that the cross section of the pair 
production of doubly charmed baryons in e + e~-annihilation could be of the same order as 
the cross section of charmonium pair production. 
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Figure 1: The leading order diagrams for the process e + e~ — > J/ip + cc. Q is momentum of the 
photon and q is momentum of the intermediate quark or gluon. The average of q 2 in the total phase 
space is about several m 2 . If both cc-pairs are fused into charmonia (e~e + — ► J/ip + charmonium), 
then q 2 ~ Q 2 /4 within the ^-approximation. 



